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Method of detecting a plurality of symbols (d k (Q) transmitted by or for a 
plurality K of users, each symbol belonging to a modulation constellation and being 
the subject of a spectral spreading by means of a spreading sequence, the said method 
comprising a filtering step (3 10,, .,,3 10*) adapted for supplying a complex vector 
(y(0.yO")) characteristic of the said received signal, the said complex vector being 
decomposed into a first vector (y*(0,y*(0) and a second vector (y'(0,y'(0 ), at least 
the closest neighbours of the first and second vectors being sought (330, 331) within a 
lattice of points (A 7 Q) generated by the said modulation constellations, the transmitted 
symbols being estimated from the components of the said closest neighbours. 



Fig 3 



The present invention concerns a multiuser detection method and device. More 
particularly, the present invention concerns a maximum likelihood multiuser detection 
method and device for a DS-CDMA (Direct Sequence Code Division Multiple Access) 
telecommunication system 
5 In a DS-CDMA mobile telecommunication system, the separation of the 

communications corning from or going to the different users is achieved by 
multiplying each complex symbol of a user by a spreading sequence which is peculiar 
to him, also referred to for this reason as the user signature. The spreading frequency 
(chip rate) being greater than the frequency of the symbols, the signal transmitted by 

10 each user is distributed (or spread) in the space of the frequencies. The ratio between 
the band occupied by the spread signal and the band occupied by the information 
signal is referred to as the spreading factor. On reception, the separation of a given 
user is obtained by means of a filtering adapted to the corresponding signature. When 
the transmission channel has a plurality of propagation paths, the output of the adapted 

15 filtering contains as many correlation peaks. Each path of the channel can be 
modelled by a complex multiplicative coefficient and a delay. The signals being 
propagated along the different paths can be aligned and combined by means of 
complex coefficients which are conjugates of the path coefficients, thus effecting a 
filtering adapted to the transmission channel. In order to simplify the terminology, the 

20 general expression "filtering adapted to the user k" will encompass both the filtering 
operation adapted to the signature of the user k and the filtering operation adapted to 
the transmission channel. 

To combat the interference between signals destined for (the downlink) or 
coming from (the uplink) the different users, multiuser detection methods have' been 

25 proposed, and notably iterative detection methods such as those known as PIC 
(Parallel Interference Cancellation) and SIC (Serial Interference Cancellation) They 
are based on the iteration of an interference elimination cycle including the estimation 
of the symbols transmitted, the evaluation of the interferences and their subtraction 
from the signals received. Although of high performance, these methods are not 

30 optimal since they do not provide an estimation in the sense of the maximum 
likelihood of the symbols transmitted by the different users. 

A method of multiuser detection with maximum likelihood inspired by the 
Viterbi algorithm was proposed by S. Verdu in an article entitled "Minimum 
probability of error for asynchronous Gaussian multiple access channels", published 



in IEEE Transactions on Information Theory, pages 85-96, January 1986, but its 
complexity is prohibitive since it varies exponentially with the number of users. 

More recently a method of multiuser detection with maximum likelihood using a 
representation by a lattice of points was proposed by L Brunei ei al , in an article 
5 entitled "Euclidian space lattice decoding for joint detection in CDMA system" 
published in Proceedings of ITW, page 129, June 1999. According to this method, a 
vector characteristic of the received signal representing a statistic sufficient for the 
maximum likelihood detection of the symbols transmitted by the different users is 
determined, it is shown under certain conditions that the characteristic vector can be 

10 represented as the point in a lattice disturbed by a noise. The detection then consists 
of seeking the point in the lattice closest to the point corresponding to the vector 
received. However, the dimension of the lattice to be used generally being 2.K or K 
where K is the number of users, the number of points to be tested is still very high. To 
simplify detection, it has been proposed to limit the search for the closest neighbour to 

1 5 the points in the lattice belonging to a sphere centred around the point received. This 
simplified detection method, referred to as the "method of detection by spheres", will 
be disclosed below: 

The context is a multiple access mobile telecommunication system with direct 
sequence spectrum spreading (DS-CDMA) comprising K users communicating 

20 synchronously with a base station. 

Let be the complex symbol sent by the user k at instant i. This symbol 
belongs to the modulation constellation A k used by the user k. which will also be 
referred to as the alphabet of symbols of the user*. 

Each user * transmits a block of N symbols with an amplitude of the signal a* 
25 The symbols are spread by a complex signature s k {t) = sf (0 +■ jsi (/) with a duration 
equal to the symbol period T : 

s k (t; = 0 if t 9. [0,T[ 

The K complex symbols U.(iydi'(i)+j.dl{i) transmitted at instant / are placed in a 

30 row vector of real values da(0 defined as. 

d:(»H^('K(^ • M^'M 0 ) 

Because of the synchronisation of the users, the corresponding modulated signal can 

be written, as a function of the time / : 
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It is assumed that the channel is an ideal channel with white additive Gaussian noise. 
Let n=S.+T}, be the signal received ai time / and 7?, a complex Gaussian noise of zeio 
mean whose components have a variance N 0 . 

Let the row vectors be y(/)=CK0,-»M')) and yii^^),yi(iX--,y»(i\yk(f)) where 
y.(i)=yk'hJy*H) ^ the complex output at instant / of the filter adapted to the user 

k: 



(3) 



— a0 

with Jtn= Ui)s t '{i)dl=Ji^J.Jtl» for M = l AT and m(()=y,..s:{t~i.T)dt 

0 0 
The autocorrelation matrix of the spreading sequences will be denoted R(i). 

If the complex elements of (3) are decomposed into their real and imaginary 
parts, there is obtained: 



1 5 Let A2=Diag(a.,a., . .,a«,a* ) and R2 be the matrix of size 2K*2K such that- 



R2 = 



-Rt.Rfr 
R&Rir 

-RM—RLRIk 



I 

(5) 
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Equation (4) can then be put in matrix form: 



y2(/)=d2(0M2-i-n2C0 



(6) 



where M2 is a real matrix of size 2K * 2K defined by M2=A 2 R 2 and where the 
noise vector M».*0M4 - ^ 0X«i(O) bas AT 0 R 2 as its covariance matrix. 



• 
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It will be demonstrated below that y 2 (0, as given by equation (6), can be 
represented as a point in a lattice A 2 of dimension 2.K, with a generator matrix IVh 

corrupted by a noise 112 

The term real lattice of points a of dimension k will be used for any set of 

vectors of R" satisfying: 

x = b,v, + b 2 Vj + ... + b^v,. where b, e Z, Vi = 1. ... ,k 
5 and where {vi,V2,. . ., var} is a base on R* 

An example of a lattice of points of dimension 2 has been shown in Fig. 1. 
The points in the lattice form an additive abelian sub-group of R* , and it is also 
the smallest sub-group of R* containing the vectors {v>,v,,...,v,} and a Z-modulus of 
R' These base vectors form the rows of the generator matrix G for the lattice. It is 
10 therefore possible to write x = bG where b={bi» . -^)eZ* . (7) 

The region delimited by the base vectors is referred to as the fundamental 
parallelotope and its volume, denoted vol(A) and det(A), is referred to as the 
fundamental volume. This fundamental volume is none other than the modulus of the 
vectorial product of the k base vectors and is therefore equal to |det(G)| where det 
15 designates the determinant. Though there are several possible choices for the 
generator matrix for the same lattice, on ihe other hand there is only one value for the 
fundamental volume. 

The Voronoi region V or Dirichlet cell of a point x belonging to the lattice is all 
the points of R* closer to x than any other point in the lattice. The volume of this 
20 region is equal to the fundamental volume. 

The stacking radius p of the lattice is the radius of the largest sphere inscribed 
in the Voronoi region and the radius of coverage that of the smallest sphere 
circumscribed in this same region. The stacking radius is therefore the radius of the 
spheres whose stacking constitutes the lattice of points and the radius of coverage is 
25 that of the smallest spheres whjch, centred on the points in the lattice, make it possible 
to cover the entire space R*. The density of the lattice is the ratio between the 
volume of the sphere of radius p and the fundamental volume. Finally, the 
coefficient of error (the kissing number) z{a) of the lattice is the number of spheres 
tangent to one and the same sphere in the stack or, in other words, the number of 

30 neighbours of a point in the lattice, situated at the minimum distance </ E = 2p 

Consider once again equation (6). The complex numbers df{i)+J.di{i) belong to a 
finite alphabet (or constellation} A of cardinal: 
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Card(A)=f\Card{Ak) (8) 

K-I 

Assume for example thai the components df{t) and are PAM modulation 
symbols of order M 

d*{i)e{-M+l-M+3r . „M-3,Af-l and (9) 
«//(/>f A/+J -Af +3,. . M-3M-1 O °) 

If the transformation is effected: 

^f(/>l(^(/>M-l) and rf , i(/>4( rf /( , ) tAf " 1 ) or a S ain vectorially; 

d' 2 (/)=l(d 2 (i)i-v^) < U > 

15 where v>* =(jVf-l ,.,,A/ -1 ) 

the components cF£ (/) and rfj(f) are elements of Z and consequently d'afc") is a vector 
ofZ* K 

In general terms, if there exists an affine transformation transforming the 
20 components d*(i) and d'(i) into elements of Z, the vector can be represented by 

a vector of Z** j 

In a similar manner, the corresponding transformation is effected on ya(/), tjiat is 

to say. 

25 y*(/)=^M'>VMM 2 ) ( 12 ) 

By means of this transformation, which is assumed to be implicit hereinafter, the 
vector d2(?)JVh then belongs to a lattice of points A2 of dimension 2.K as defined by 
equation (7) with G=Mi . The vector y 2 (0 can then be considered 10 be a point in the 
30 lattice A 2 corrupted by a noise n 2 (i). 

If 11 is assumed that the components of the noise vector n 2 (0 are centred 
independent random Gaussian variables, the problem of detection in the sense of the 
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maximum likelihood of the symbols transmitted by the different users amounts to a 
search for the point z 2 in the lattice A 2 such that its distance to y 2 (0 is at a minimum. 

In reality, the components of the noise vector n 2 (0 are correlated and the 
covariance matrix of n 2 (0 is N 0 R 2 . 
5 In order 10 boil down to the decollated case n is necessary, prior to the 

decoding, to effect an operation of whitening of the noise. 

The matrix R being hermitian, the autocorrelation matrix R 2 is symmetrical 
defined positive and can therefore be the subject of a Cholesky factorisation: 

R 2 =W 2 W/ < 13 > 
] 0 where W? is an inferior triangular matrix of size 2Kx2X. 

A whitened observation vector is defined: yt{i)*y2{i}Wf~ l (14) 
and a new lattice of points fi 2 consisting of vectors of components 
-MfcM wilh Hh^M' 1 vvhere x 2 (i) is a vector of components 

(*i*WW(i • .jftixilt)) belonging to A 2 . 

j 5 i t can easily be shown that, after whitening, the covariance matrix of the filtered 

noise nsOJWf"' is equal to Milz* where 1 2 k is the identity matrix of dimension 2K. 
The detection therefore comprises a first step of whitening the observation vector 
followed by a step of seeking the closest neighbour within the lattice of points Q 2 . 

In order to reduce the number of points to be tested, as illustrated in Fig. 1 , it is 

20 possible to limit the search to a sphere centred around the point y 2 . In practice, the 
choice of the radius of the sphere results from a compromise: n must not be too large 
in order not to lead to an excessively high number of points and sufficiently large to 
include at least the closest neighbour. 

Fig. 2 depicts schematically a multiuser detection device using a method of 

25 detection by spheres. The received signal n is filtered by a battery of filters adapted 
to each of the users, 210 I ,... T 210 K . The real and imaginary components of the 
observation vector y 2 (0 output from the adapted filters are transmitted to a matrix 
calculation unit performing the spectral whitening operation according to equation 
(14). The real and imaginary components of the whitened vector y2(i) are then 

30 transmitted to a unit for detection by spheres seeking the closest neighbour of the 
point received within the lattice Q 2 of dimension 2./C The coordinates of the closest 
neighbour directly give the real and imaginary components of the estimated symbols 
dk{i) for the different users 
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It can be shown that the method of detection by spheres as disclosed above has a 
complexity in terms of 0(K") which may prove very disadvantageous when the 

number of users is high 

The aim of the present invention is to propose, under certain conditions, a 
5 simplification of the method of detection by spheres. 

For this purpose, the invention is defined by a method of detecting a plurality of 
symbols (d k (/)) transmitted by or for a plurality K of users, each symbol belonging to a 
modulation constellation and being the subject of a spectral spreading by a spreading 
sequence, the said method comprising a filtering step adapted for supplying a complex 
10 vector (y(O,y(0) characteristic of the said received signal, the said complex vector 
being decomposed into a first vector (y\i)S"(0) * second vector (y'(0,?'(0) 
and at least the closest neighbours of the first and second vectors being sought within 
a lattice of points (A.Q) generated by the said modulation constellations, the 
transmitted symbols being estimated from the components of the said closest 

15 neighbours o 

Advantageously, the spreading sequences ( s t {t) ) consist of real multiples ( *(0 ) 

of the same complex coefficient (a). 

According to one embodiment of the invention, the search is limited to a first set 
of points in the lattice belonging to a first predetermined zone (Ir) around the first 
20 vector and a second set of points in the lattice belonging to a second predetermined 
zone (Zj) around the second vector. 

Likewise ihe search can be limited to a first set of points in the lattice belonging 
to a first predetermined zone (£r) around the origin and a second set of points m the 
lattice belonging to a second predetermined zone (Z0 around the origin. I 
25 The said first and second predetermined zones are for example spheres. 

Advantageously, the search for the closest neighbour of the first vector is carried 
out on a plurality of components of the latter, the search being limited for each of the 
said components to an interval defined by a lower band and an upper band, the said 
bands being chosen so that the said interval does not comprise any points relating to 
30 symbols which cannot belong to the modulation constellation. 

Likewise the search for the closest neighbour of the second vector can be carried 
out on a plurality of components of the latter, the search being limited for each of the 
said components to an interval defined by a lower bound and an upper bound, the said 
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bounds being chosen so that the said interval does not comprise any points relating to 
symbols which cannot belong to the modulation constellation. 

Prior to the search for the closest neighbour, the first vector (y*(i)) »s 
advantageously subjected to a matrix processing aimed at substantially decollating 
the different noise components thereof 

Likewise, prior to the search for the closest neighbour, the second vector ( y'(0) 
can be subjected to a matrix processing aimed at substantially decorrelating the 
different noise components thereof. 

According to a variant embodiment of the invention, the said search step is 
extended to a search for a first set of points which are closest neighbours of the said 
first vector, referred to as first neighbours, and to a second set of pomts which are 
closest neighbours of the said second vector, referred to as second neighbours, and the 
transmitted symbols are estimated flexibly from the symbols generating the said first 
and second neighbours and distances separating the said first neighbours of said first 
vector on the one hand and the said second neighbours of the said second vector on 
the other hand. 

According to a particular embodiment of the invention, there are determined, 
from the estimated symbols, the contributions of each user to the signals obtained by 
the adapted filtering step and, for a given user k 9 the contributions of the other users 
corresponding to the symbols already estimated are eliminated at the output of the 
filtering step. Alternatively, the contributions of each user to the signal received from 
ihe estimated symbols are determined and, for a given user k, the contributions of the 
other users corresponding to the symbols already estimated are eliminated at the input 
of the adapted filtering stage. 

If the symbols of the said K users are transmitted synchronously, the said lattice 
of points will preferably be of dimension K, 

If the symbols of the said K users are transmitted asynchronously and propagate 
along a plurality of paths, the dimension of the lattice will preferably be equal to the 
number of symbols of the different users which may interfere and are not yet 
estimated. 

The present invention is also defined by a device for detecting a plurality of 
symbols (d k (0) transmitted by or for a plurality K of users, each symbol belonging to a 
modulation constellation and being the subjeci of a spectral spreading by means of a 
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spreading sequence, the device comprising means for implementing the method 
disclosed above. 

Tnis device can in particular be used in a receiver in a DS-CDMA mobile 
lelecommunjcations system. 
5 The characteristics of the invention mentioned above, as well as others, will 

emerge more clearly from a reading of the following description given in relation to 
the accompanying figures, amongst which: 

Fig. 1 depicts a lattice of points useful to the detection method employed in the 
receiver illustrated in Fig. 2; 
l0 pig 2 depicts schematically the structure of a multiuser DS-CDMA receiver 

using a method of detection by spheres; 

Fig. 3 depicts schematically the structure of a multiuser detection device 
according to a first embodiment of the invention; 

Fig. 4 depicts schematically the structure of a multiuser detection device 
15 according to a second embodiment of the invention. 

The idea at the basis of the invention is to reduce the dimension of the lattice of 
points using spreading sequences of a particular type for the different users. 

Consider once again a DS-CDMA telecommunication system with K 
synchronous users. If spreading sequences *(r) with real values are chosen, the 
20 imaginary terms of the matrix R 2 and consequently of the matrix Mitf are zero. 
Consequently, the system can be modelled by a lattice of real points A of dimension 
K and generator matrix M(j) : j 

f 

v*(z)=d*(0M(i>n*W i l 5 ) 

25 y'(0=d'(OM(0-i-n'(7) < l6) 

where y*(0,d*(/),n*« (or respectively y'(0,d'(0,n'(O) ar e the vectors consisting of 
the real pans (or respectively of the imaginary pans) of the components of 
y(/),d(/),n(/) ; 

T 

30 M(/)=AR(0 where R(i) is the matrix consisting of the coefficients A«« pf(r)si(0dr 

u 

and A is the vector of the amplitudes of the K users. 
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The observation vectors y*(/) and y'(0 belong to R x . After any transformation 
according to an equation of the same type as (12), the vectors y*(i) and y'(i) can be 
considered to be points in a lattice A of generator matrix M(0 corrupted by noise 

It can easily be shown that the noise vectors n*(0 and n'(/) both have the 
covariance matrix No.R(0 » bein 8 a symmetrical matrix defined positive, it can be 
facionsed according to a Cholesky decomposition: R «= WW T where W is an 
inferior triangular real matrix of size KxK In order to decorrelate the noise 
components, the real observation vectors y*(0 and y'(0 are first of all subjected to a 
whitening operation: 

y'OHW' (17) 

nwdjwT (18) 

Secondly, the closest neighbours of the vectors y*(0 and y'(j) belonging to the 
1 5 lattice of points Cl consisting of the vectors xO>x(OW T ~' where x(i) belongs to A, are 
sought. It can easily be shown that, after whitening, the covariance matrices of the 
filtered noises o'OJw 7 "' and n'OJW 7- ' are both equal to MI* where I* is the 

identity matrix of dimension K. 

It can therefore be seen that using real signatures leads to a search for two closer 
20 neighbours in the same lattice of dimension K whilst in the general case, which is 
complex, the decoding requires a search in a lattice of dimension 2K. In fact it can 
easily be shown that the above result extends to any set of signatures carried by the 
same complex number, that is to say such that: where a is a complex 

number and s°(t) is real. This will of course in particular be the case if all the 

25 signatures are imaginary. 

Since the search for the closest neighbour for the vector y'(;) takes place 
according to the same principle as for y*(z), only the first will be disclosed This 
search consists of determining the point x minimising the metric: 

4>7x}=iKH ? +Hf (19) 

30 where y"=x+n" and x=(jr., . ,xr) is a point belonging to the lattice fi. 

Alternatively, it will be noted that the vector y\i) does not need to be whitened 
if a metnc is used based on the covariance matrix 
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mi?Vx}*{y*-x)R-Xy*-x)' (20) 

Hereinafter, for reasons of simplification, the observation vector, whitened 
5 ( y A (0 ) or not ( y*(0 ), will be termed z and the metric acting in equation (1 9) or (20) 
will be termed |j| • 

The points in the lattice Ci consist of the vectors x such that x=bG where G is 
the generator matrix for the lattice and b=(h,...,&), the components h belong to the 

ring of integers Z . 

!0 The detector advantageously restricts its metric calculation to the points which 

are situated within a zone of the constellation situated around the received point, 
preferably within a sphere of given radius Vc centred on the received point z. Only 
the points in the lattice situated at a quadratic distance less than C from the received 
point are therefore considered for the minimisation of die metric. 

15 In practice, the decoder effects the following minimisation: 

"S^isSM (21) 

To do this, the decoder seeks the smallest vector w in the translated set z-ft . 
20 The vectors z and w can be expressed as: 

z^pG with p=(/a>-->P0 j 
w=^G with f 2) 

25 It is important to note that p and 4 are real vectors. As w=*z-x, where x 

belongs to the lattice Q, this gives the equation £ =p t -b t for /=l r .. 7 *r with 

The vector w is a point in the lattice whose coordinates £ are expressed in the 
translated reference frame centred on the received point. The vector w belongs to a 
30 sphere of quadratic radius C centred at 0 if. 



(23) 
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In the new system of coordinates defined by £, . the sphere of quadratic radhis C 
centred ai y is therefore transformed into an ellipsoid centred on the origin. The 
Cholesky factorisation of the Gram matrix V=GC T gives r=AA r , where A is an 
inferior triangular matrix of elements S(, . 

It should be noted that, if the vector y has been whitened, it is not necessary to 
effect this factorisation since the generator matrix of CI is equal to AW and is 
therefore already triangular and inferior. However, where prior whitening has not 
been carried out, Cholesky decomposition is necessary. In all cases, it is possible to 
write: 

otO^V+^l 1 -^***^**!^ (24) 

By putting 

q u m S.. 2 for 1 = 1, ... ,K 

q.. = ^- for j = 1, ... ,k; i = j f l, ••• ,k 

^ jj 

15 there is obtained 

qgo-&4**£«*] (25) 

By taking firsi of all the range of possible variations of & , and then adding the 
components one by one, the following K inequalities are obtained, which define all the 
20 points within the ellipse: 

q,-i.*-^+q*.*-& f+q^ZC ( 26 ) 
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It can be shown that the inequalities (26) make it necessary for the integer components 
of b to satisfy: 
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+ A + ]£ff>ft 



(27) 



Si. 



C - X 9« ft + £ */rft + A + Z <lj£j 



where |Y| is the smallest integer greater than the real number x and [xj is the largest 

5 integer smaller than the real number x . 

The decoder has K internal counters, namely one counter per dimension, each 
counter counting between a lower and upper bound as indicated in (27), given that 
each counter is associated with a particular pair of bounds. In practice these bounds 
can be updated recursively. 
10 Advantageously, all the values of the vector b are listed for which the 

corresponding point in the lattice x = bG is situated within the quadratic distance C 
from the received point. The points in the lattice situated outside the sphere in 
question are not tested. Advantageously, for each component the upper and 

lower bounds of the search will be adjusted so as not to contain points which are 
certainly outside the constellation. Thus the counters do not waste any time in 
running over points which, in any event, are not solutions. For example, if all the 
users employ the same PAM modulation constellation of si2e A/ T the search bounds 
cannot leave the interval [0,M-1]. 

In addition, the search within the sphere can be considerably accelerated by 
20 updating the radius with the last Euclidian norm calculated ||w||. Finally, there is 
selected, as the best point x , the one associated with the smallest norm ||w|| . 

So as to be sure that the decoder finds at least one point in the lattice, a search 
radius is advantageously chosen so as to be greater than the radius of coverage of the 
lattice. It can for example be taken so as to be equal to the upper Rogers bound: 
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given instant may interfere with two successive spread symbols of another user k\ If 
it is assumed that the scattering of the transmission delays t a of the different users is 
less than a symbol period T t the symbol of a user k sem at an instant t 9 d k (i), may 
interfere with the symbols d k {Ul) and d k (i+l) of a user k\ It can be assumed without 
loss of generality that 0<r»< . <r*< . <tk<T . After the symbols d^i-1) for all the users 
k=J~ K are detected, the detection is commenced of the symbols relating 10 the instant 
i commencing with the earliest user (here the user 1) and finishing with the latest user 
(here the user K). The detection of d k (i) depends on three vectors- the vector 
dp=(rfi(O f ...//i-i(0,tf*('- 1 )» A '»('-l)»- relating to the symbols already 
detected, the vector dF^i/i(i>l),...^«H(f>l)^(0^'ri(Oi--/'^(0) relating 10 future 
symbols and the observation vector y(0 . It can be shown that the complex 
observation vector can be written in the form of a past contribution and a future 
contribution; 

y=drARpi-d fARf+d (30) 
where n is a noise vector of covariance matrix A^Rp, and Rp and Rp are respectively 
the matrix of correlation of the signatures with the passed signatures and the matrix of 
correlation of the signatures with the future signatures. 

When the signatures are of the real value type (or more generally are real 
multiples of the same complex number) equation (30) can be decomposed in the form 
of two equations involving only real vectors, in a similar manner to equations (15) and 
(16): 

y*=dF*ARp+dp*ARF+n* (31) 
y^p'ARpWARp-hn' (32) 

Instead of effecting the detection by spheres on the observation vectors y* and 
y\ this will operate on the vectors freed of any interference due to the passed symbols 
already estimated, namely: y*-dp*ARp and y'-dr'ARr. 

Alternatively, instead of carrying out the subtractive elimination of the 
interference at the level of the observation vector, this elimination can be envisaged in 
an equivalent manner, upstream, on the spread of signals- The estimated symbols are 
respread spectrally and the contributions of the different users to the received signal 
are subtracted one by one, on the fly (as they are respread), at the input of the adapted 
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VC^(/Clo^+/Clog]ogA:+5A:>^^ (28) 

where Vk is the volume of a sphere of unity radius in rhe real space. 

Fig. 3 illustrates schematically a multiuser detection device according to a first 
embodiment of the invention The received signal is first of all filtered by a battery of 
fillers adapted to the different users, 310!,. ,310*. The observation vector output from 
the adapted filters is decomposed into a real observation vector y*(0K^*(0» • >>*(0) 
and an imaginary observation vector y'(0K>A0>- J^CO) • After any transformation of 
the type at (12) (not shown), the vectors y"(i) and y'(z') undergo a spectral whitening 
at 320 and 321 in order to decorrelate the noise samples. The whitened vectors y*(/) 
and y'(i) are then the subject of a search for the closest neighbour as described above 
in the detectors by spheres 330 and 331. The point found by the detector 330 gives 
(by means if necessary of a transformation which is the reverse of that of (12)) the real 
components of the symbols estimated for the K users. Likewise, the detector 331 
gives the imaginary components of these estimated symbols. 

Instead of directly supplying symbols of the constellation, the receiver can be 
adapted to supply symbols in the form of flexible decisions. In this case, the search 
within the detection sphere is no longer limited to the closest neighbour but is 
extended to a plurality of the closest neighbours of the point relating to the received 
signal. 

More precisely, let 2* and 2, be the spheres centred respectively around y*(i) 
and y'(/). There is associated, with any pair d m "(OK of neighbouring points 

belonging to (2 Rl I,) such that the K components d? M [i) belong to the modulation 
constellations of the users, an a posteriori probability p'"' M , a probability that the 
vector rff '"(0 defined by this point has been sent, given the observation %), Let © 
be the set of these pairs. A flexible symbol of a user k is defined as the Mriuplet 
(*/,.... Xmk) where M k is the cardinal of the modulation constellation of the user k and 
where 75 is the probability that the symbol s, has been sent. This gives: 

2>^ /rf "^^ (29) 
The a posteriori probabilities p m m can for example be expressed as a function of 

the distances A„ t and X m - separating the vectors y*(i) and y'{i) from v m and v m . 

Equation (2) presupposed that the signals of the different users were 

synchronous. When this assumption is not valid, the spread symbol of a user k at a 
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filters The inputs of the adapted fillers will therefore receive the signals cleaned of 
the passed contributions, namely n-h(t) where: 

/aCO^^'-COM'-^ (33) 

5 

The first right-hand term in equation (33) represents the contribution of the users 
which are in advance on the user k and for which consequently an estimation of the 
symbols dr(t) is already available. The second term represents the contributions of 
the users which are behind with respect to the user k: the contribution due to the 

10 symbols dk(i) can of course not be evaluated at the time when the symbol dt(t) is 
despread. On the other hand, the contribution due to the previous symbols A</-1) 
can be evaluated since an estimation of the symbols J*</-I) is already available. 

It is important to note that, unlike the synchronous case where the step of 
detection by spheres estimated only all the symbols dk(i) , in the asynchronous case 

15 the symbols of the different users are estimated one after the other in their order of 
arrival 

Fig. 4 illustrates a second embodiment of the invention implementing a 
subiractive elimination of the interference for asynchronous users. To simplify the 
representation, only the detection branch concerning a user k has been illustrated. The 

20 sum symbolises the sum of the contributions of the other users and is subtracted 
at 405* from the signal n input from the adapted filter 410,. The detection by spheres 
for the user k starts at the instant when the symbol dk(J) is despread. It should be 
noted that, at this instant, the interfering symbols not yet estimated of the other users 
will have been subject only to a partial despreading (except of course if the interfering 

25 symbols are synchronous. The module 430* (or respectively 431 a) operates on the 
components output from 42<\ (or respectively 42 U). Alternatively the outputs of all 
the modules 430* (or respectively 43 1 A ) can be demultiplexed in order to be processed 
by a common module 440 working K times more rapidly. Unlike Fig. 3, here only the 
ouLput of the modules 430* and 43 U is used for supplying the estimated symbol 

30 d>(.<) . The symbol Mi) is then respread and the contribution akdi(i)si(t-iT-n) of 
ihe user k is evaluated and then subtracted at the input of the adapted filters 410 A - via 

If the transmission channels of the different users are of the multipath type, the 
problem of the elimination of interference is more complex since it is necessary to 
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consider the interference between all the paths of all the users. The received signal 
can be written 

" = S£]>- C * <>< **(' M'-iT-Ti* h-QAO (34) 

,-0 km] pm\ 

where P k is the number of parts of the transmission channel of the user k M Tp*=Ti+0j* 
is the total delay of the delay (r* ) on transmission of the user k and the delay (ft* ) of 
propagation along the path p of the channel k and c p a is the complex multiplicative 
factor associated with this path. It 'will be assumed once again that 
0<tl<..<ta< , <t*<T and that in addition the scattering of the paths is less than the 

symbol period: 0<>6p*<T . As a result Q<n<2T . 

The adapted filters 41 0 k then effect a filtering adapted to the signature of the user 

k and to the transmission channel k by MRC {Maximum Ratio Combining) of the 

signals of the different paths. More precisely, the filter 310 k performs the following 

operation; 

M0=2_^>* \sk\i-iT~TpL)nto (35) 

Because of the scattering of the delays, the duration to be taken into account for 
the interference between users is 2T and the symbol of a user can interfere with two 
consecutive symbols, not yet estimated, of another user. The conjoint detection then 
relates to all the interfering symbols which have not yet been estimated, the estimated 
interfering symbols for their part being used for the subtractive elimination. The 
complexity of the detector by spheres is higher than in the single-path asynchronous 
case since the dimension of the lattice to be taken into account is 2/f-l instead of AT. 

Finally, it should also be noted that the estimation of the interfering 
contributions in the single-path or multi-path asynchronous case may be effected 
either from constellation symbols (hard decisions) or from the flexible symbols of the 
different users. 

Although certain embodiments of the invention have been depicted in the form 
of functional modules, it is clear that the device according to the invention can be 
produced in the form of a processor programmed for executing the different functions 
depicted or in the form of a plurality of dedicated processors able to implement one or 
more of these functions. 
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CLAIMS 

1. Method of detecting a plurality of symbols (d k (/)) transmitted by or for a 
plurality K of users, each symbol belonging to a modulation constellation and being 
the subject of a spectral spreading by means of a spreading sequence, the said method 
comprising a filtering step (310,,..,310 A ) adapted for supplying a complex vector 

5 (y(0,y(0) characteristic of the said received signal, characterised in that the said 
complex vector is decomposed into a first vector ( y*(0,y*CO ) and a' second vector 
(y'(0,y'(/)) and in that at least the closest neighbours of the first and second vectors 
are sought (330, 331) within a lattice of points (A,Q) generated by the said modulation 
constellations, the transmitted symbols being estimated from the components of the 
1 0 said closest neighbours. 

2. Detection method according to Claim 1, characterised in that the spreading 
sequences {s A (t)) consist of real multiples (*?(')) of the same complex coefficient 

15 

3. Detection method according to Claim 1 or 2, characterised in thai the search 
is limited to a first set of points in the lattice belonging to a first predetermined zone 
(£ R ) around the first vector and a second set of points in the lattice belonging to a 
second predetermined zone (Ii) around the second vector. 

20 J 

4. Detection method according to Claim 1 or 2, characterised in that the search 
is limned to a first set of points in the lattice belonging to a first predetermine^ zone 

around the origin and a second set of points in the lattice belonging to a second 
predetermined 2one (Z,) around the origin. 

25 

5. Detection method according to Claim 3 or 4, characterised in that the said 
first and second predetermined zones are spheres. 

6 Detection method according to one of the preceding claims, characterised m 
30 that the search for the closest neighbour of the first vector is effected on a plurality of 
components thereof, the search being limited for each of the said components to an 
interval defined for a lower bound and an upper bound, the said bounds being chosen 
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so thai the said interval does not comprise points relating to symbols which cannot 
belong to the modulation constellation. 

7. Detection method according to one of the preceding claims, characterised in 
that the search for the closest neighbour of the second vector is effected on a plurality 
of components thereof, the search being limited for each of the said components to 
an interval defined for a lower bound and an upper bound, the said bounds being 
chosen so that the said interval does not comprise points relating to symbols which 
cannot belong to the modulation constellation. 

8. Detection method according to one of the preceding claims, characterised in 
that, prior to the search for the closest neighbour, the first vector ( y*(0 ) is subjected 
to a matrix processing (320) aimed at substantially decollating the different noise 
components thereof. 

9. Detection method according to one of the preceding claims, characterised in 
that, prior to the search for the closest neighbour, the second vector (y'(/)) is 
subjected to a matrix processing (321) aimed at substantially decoirelating the 
different noise components thereof. 



10. Detection method according to one of the preceding claims, characterised 
in that the said search step is extended to the search for a first set of points which are 
the closest neighbours of the said first vector, referred to as first neighbours, and a 
second set of points which are closest to the said second vector, referred to as second 
25 neighbours, and in that the transmitted symbols are estimated flexibly from symbols 
generating the said first and second neighbours and distances separating the said first 
neighbours from the first vector on the one hand and the said second neighbours 
from the said second vector on the other hand. 

30 11. Detection method according to one of the preceding claims, characterised 

in thai the contributions of each user to the signals obtained by the adapted filtering 
step are determined from the estimated symbols and in that, for a given user k, the 
contributions of the other users corresponding to the symbols already estimated are 
eliminated at the output of the filtering step. 
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12. Detection method according to one of Claims 1 10 10, characterised in that 
ihe contributions of each user to the received signal are determined (340) from the 
estimated symbols and in that, for a given user k, the contributions of the other users 

5 corresponding to the symbols already estimated are eliminated at the input of the 
filtering step. 

13. Detection method according to one of Claims 1 to 10, characterised in that, 
ihe symbols of the said K users being transmitted synchronously, the said lattice of 

10 points is of dimension K. 

14. Detection method according to Claim 11 or 12, characterised m that, the 
symbols of the said K users being transmitted asynchronously and propagating along 
a plurality of paths, the dimension of the lattice is equal to the number of symbols of 

15 the different users which may interfere and are not yet estimated. 

15. Device for detecting a plurality of symbols (d k (i)) transmitted by or for a 
plurality K of users, each symbol belonging to a modulation constellation and being 
the subject of a spectral spreading by a spreading sequence, the device comprising 

20 means for implementing the method claimed according to one of the preceding 
claims. 

I 

16. Receiver for a DS-CDMA mobile telecommunication system comprising a 
detection device according to Claim 15. \ 
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ABSTRACT 

Multiuser detection method and device in DS-CDMA mode 

Method of detecting a plurality of symbols (d k (0) transmitted by or for a 
plurality K of users, each symbol belonging to a modulation constellation and being 
the subject of a spectral spreading by means of a spreading sequence, the said method 
comprising a filtering step (310,, ..,310*) adapted for supplying a complex vector 
(y(i),y(0) characteristic of the said received signal, the said complex vector being 
decomposed into a first vector (y*(/),y*(0 ) and a second vector (y ; (0»y'(0 \ at least 
the closest neighbours of the first and second vectors being sought (330, 331) within a 
lattice of points (A,QJ generated by the said modulation constellations, the transmitted 
symbols being estimated from the components of The said closest neighbours. 

Fig. 3 



